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Abstract 

Core collapse supernovae are a huge source of all flavor neutrinos. The flavor composition, 
energy spectrum and time structure of the neutrino burst from a galactic supernova can 
provide information about the explosion mechanism and the mechanisms of proto neutron 
star cooling. Such data can also give information about the intrinsic properties of the 
neutrino such as flavor oscillations. One important question is to understand to which 
extend can the supernova and the neutrino physics be decoupled in the observation of a 
single supernova. On one hand, the understanding of the supernova explosion mechanism 
is still plagued by uncertainties which have an impact on the precision with which one 
can predict time, energy and flavor-dependent neutrino fluxes. On the other hand, the 
neutrino mixing properties are not fully known, since the type of mass hierarchy and 
the value of the ^13 angle are unknown, and in fact large uncertainty still exists on 
the prediction of the actual effect of neutrino oscillations in the event of a supernova 
explosion. In this paper we discuss the possibility to probe the neutrino mixing angle 613 
and the type of mass hierarchy from the detection of supernova neutrinos with a liquid 
argon TPC detector. Moreover, describing the supernova neutrino emission by a set of 
five parameters (average energy of the different neutrino flavors, their relative luminosity 
and the total supernova binding energy), we quantitatively study how it is possible to 
constrain these parameters. A characteristic feature of the liquid argon TPC is the 
accessibility to four independent detection channels ( (1) elastic scattering off electrons, 
(2) charged neutrino and (3) antineutrino and (4) neutral currents on argon nuclei) which 
have different sensitivities to electron-neutrino, anti-electron-neutrino and other neutrino 
flavors (muon and tau (anti)neutrinos). This allows to over-constrain the five supernova 
and the fiavor mixing parameters and to some extent disentangle neutrino from supernova 
physics. Numerically, we find that a very massive liquid argon detector (0(IOO kton)) 
is needed to perform accurate measurements of these parameters, specially in supernova 
scenarios where the average energies of electron and non-electron neutrinos are similar 
(almost degenerate neutrinos) or if no information about the ^13 mixing angle and type 
of mass hierarchy is available. 
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1 Introduction 



Core collapse supernovac arc a huge source of all flavor neutrinos. Neutrino astrophysics en- 
tered a new phase with the detection of neutrinos from the supernova SN1987A in the Large 
Magellanic Cloud by the Kamiokande and 1MB detectors [1]. In spite of these fundamental 
neutrino observations, the 19 events observed are not statistically significant enough to ob- 
tain precise quantitative information on the neutrino spectrum. Currently running neutrino 
detectors like Supcrkamiokandc [2] or SNO [3] have the capabilities to provide high statistics 
information about supernova and neutrino properties if a supernova collapse were to take 
place in the near future. 

The flavor composition, energy spectrum and time structure of the neutrino burst from a 
galactic supernova can give information about the explosion mechanism and the mechanisms 
of proto neutron star cooling. 

The neutrino signal from a galactic supernova can also give information about the flavor 
oscillation of neutrinos. Although new data from solar, atmospheric, reactor and accelerator 
neutrinos [4] have contributed to the understanding of the neutrino properties, still the neu- 
trino mixing angle ^13 and the nature of the mass hierarchy, i.e., normal or inverted, remain 
unknown. Only an upper limit on the mixing angle is available dominated by negative mea- 
surements from reactor experiments [5]. These parameters can be probed by the observation 
of supernovae neutrino bursts, since neutrinos will travel long distances before reaching the 
Earth and will, as they travel through the exploding star, in addition traverse regions of 
different matter densities where matter enhanced oscillations [6] will take place. In the case 
of the angle, matter enhancement has the striking feature that very small mixing angles, 
beyond any value detectable by next generation accelerators, could in fact alter significantly 
the ncTitrino spectrum. Hence, supernova neutrinos provide indeed a complementary tool to 
study the 6*13 angle. 

The main question of course is to understand to which extend can the supernova and the 
neutrino physics be decoupled in the observation of a single supernova. On one hand, the 
understanding of the supernova explosion mechanism is still plagued by uncertainties which 
have an impact on the precision with which one can predict time, energy and flavor-dependent 
neutrino fluxes. On the other hand, the intrinsic neutrino properties are not fully known, 
since the type of mass hierarchy and the value of the ^13 angle are unknown. 

The total event rates and the energy spectra of the detected supernova neutrinos depend 
on the neutrino flavor and on the ^13 mixing angle and the type of mass hierarchy. However, 
it is not a priori easy to separate neutrino and supernova physics. Different solutions to this 
problem have been proposed. Global fits of data ([7, 8]) and analyses considering qualitative 
features of the supernova neutrino fluxes ([9]) have been performed for the case of Cerenkov 
detectors. The potential of liquid scintillator detectors has also be investigated ([10]). 

We have already qualitatively illustrated in our paper [11] that four independent channels 
are a priori accessible with a liquid argon TPC: elastic scattering on electrons, charged 
current absorption on argon, charged current absorption on argon and neutral current 
interactions on argon. They have different sensitivities to electron-neutrino, anti-electron- 
neutrino and other neutrino flavors (muon and tau (anti)neutrinos). These events being 
sensitive to different combinations of neutrino and antineutrino flavors can be combined to 
over-constrain the features of the supernova fluxes that are at the moment not well constrained 
by theory and the mixing parameters and allow neutrino and supernova physics effects to be 
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disentangled. 

In the present paper we study quantitatively the capability of LAr TPC detectors to 
simultaneously measure the supernova and oscillation parameters, performing a global fit to 
the event energy distribution for elastic and charged current events and to the rate for the 
neutral current events. 

This study is motivated by the 3 kton LAr ICARUS experiment [12] that was proposed 
at the Gran Sasso laboratory, and by ideas for a very large liquid argon TPC with a mass in 
the range of 100 ktons (see [13] and independently [14]). 

2 Theoretical framework 

We briefly summarize in this section the neutrino features relevant for this analysis. We 
essentially follow the general framework already described in our previous paper [11]. 

2.1 Supernova physics 

At the end of the life of a massive star with a core mass greater than the Chandrasekhar 
mass, the energy generation in the core stops, and the star undergoes a fast gravitational 
collapse. Within a duration of tens of milliseconds, this collapse compresses the inner core 
beyond nuclear matter densities. During this stage electron neutrinos are produced via neu- 
tronisation processes. While the electron neutrinos produced could escape the star during the 
first period (producing the "neutrino burst"), the very dense core becomes opaque to neu- 
trinos which remain captured within the supernova core producing the so-called inner core 
rebound. This generates a shock wave which propagates radially outward through the star 
remnants, supplying the "explosion" mechanism. During the shock phase, thermal neutrinos 
and antineutrinos of all flavors are produced due to high temperatures. It is believed that 
the 99% of the total binding energy of the star, Eb ~ 3 x 10^^ ergs, is emitted in the form 
of neutrinos. These neutrinos are in equilibrium with their surrounding matter density and 
their energy spectra can be described by a function close to a Fermi-Dirac distribution. The 
flux of a neutrino Va emitted can then be written as [9] : 

La 

where is the luminosity of the flavor {Eb= Y1 ^a)-, D is the distance to the supernova, 
Eoi is the energy of the neutrino, is the neutrino temperature inside the neutrinosphere 
and ?7a is the "pinching" factor. We have taken ??a = and hence no deviations in the high 
energy tail are considered. The neutrino average energy {Ey^) depends on both Tq, and ry^. 
For -qa = 0, {Ey^) k, 3.15 T^. The normalization factor in absence of pinching is equal to 
F3(0) ^ 5.68. 

The original v^, u^-, 9^ and i>T fluxes are approximately equal and therefore we treat them 
as Vx- Out of all neutrinos, muon and tau neutrinos and their antiparticles interact with 
matter through neutral currents only, whereas electron neutrinos and their antiparticles feel 
both charged and neutral currents. Therefore, muon and tau flavors should decouple first and 
hence have the largest temperatures. This energy hierarchy between the different neutrino 
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flavors is generally believed to hold and imply {E^^) < (Ep^) < (E,,^), due to the different 
strength of neutrino interactions with the surrounding matter. However, the specific neutrino 
spectra remain a matter of detailed calculations. In particular, the simulations of Raffelt et 
al. [15] seem to indicate that the energy differences between flavors could be very small. 
We use a range of typical values of the neutrino average energies and relative luminosities 
provided by several simulations in the following intervals [9, 15]: 

{E^J= (7-18) MeV; {EpJ= (14-22) MeV; {E^^)= (15-35) MeV (2) 
Le/L^= (0.5-2); = Lg (3) 

In practice, we will consider two speciflc scenarios (I &: II) in order to understand the effects 
of a hierarchical versus non-hierarchical distribution of energies on our results, see below, 
and take luminosity equipartition. Some detailed simulations may also give rise to cases 
where the electron luminosities of the neutrinos vary between (0.5-2) times the muon and 
tau luminosities. We did not explicitly consider this case, even though our calculations could 
be extended to such a case. We will be later concerned on how well one can determine 
the energies and the relative luminosities of the neutrinos through the detection of a single 
supernova collapse. 

We consider core collapse supernovae at a distance D of 10 kpc (galactic supernovae). 

We assume that in the actual event of a supernova, this distance will be precisely estimated 
by independent astronomical observations and we hence take this value as precisely known 
in our discussion^. 

In order to parameterize the supernova physics, we employ five supernova parameters: 

1. the total binding energy Eb= Y1 j 

2. the average energies of the neutrinos emitted from the supernova {Ei,J, {EpJ, {Ey^) and 

3. the relative luminosities of the electron and non-electron neutrinos L^/Lx- 

We always assume Lg = Lg and = L^, hence we obtain: 

Eb = Lc + Le + 2Lx + 2Lx = 2Le + AL^ 

= 2Lx{LjLx + 2) (4) 

from which = = Eb / {2{Le/ + 2)). 

In order to take into account the uncertainties on the knowledge of the supernova physics, 
we consider two different sets of reference values for the supernova parameters from Ref. [16] 
and [15], corresponding to one hierarchical energy scenario (I) and a non-hierarchical energy 
scenario (II). They differ in the assumptions on the average energies of the neutrino flavors 
(Table 1). 

^We note that the error on the distance will essentially affect the overall normalization of the fluxes, 
parameterized by the binding energy iu our discussion. This error will not fundamentally affect our conclusions. 

On the other hand, the fact that a suijernova could actually occur at a closer or larger distance does affect 
the statistics of collected event. We take 10 kpc as a reference value for easy comparison with other authors. 
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SN 


Eb 








Luminosity 


Ref. 


scenario 


(x 10^3 erg) 


(MeV) 


(MeV) 


(MeV) 






I 


3 


11 


16 


25 




[16] 


II 


3 


13 


16 


17.6 




[15] 



Table 1: Assumed supernova parameters: scenarios I and II. 



2.2 Neutrino flavor oscillation physics 

We summarize here the supernova flux changes introduced by flavor oscillations (sec e.g. 
Ref. [11] for more details). Before arriving at Earth the neutrino fluxes travel through the 
supernova mantle undergoing two MSW [6] resonances: one at high density p k, 10^-10^ g 
cm~^ (H-resonance) which is governed by the atmospheric parameters {^m\i and ^13) and 
the other one at low density p ~ 10-30 g cm~^ (L-resonance) , characterized by the solar 
parameters (Am2i and 9i2)- The exact propagation of the neutrinos and their oscillation 
from one flavor to another depend on the matter density profile. We refer the reader to Ref. 
[11] for more details on our implementation of these effects. 

The transitions in the two resonance layers can be considered independently and each 
transition is reduced to a two neutrino oscillation problem. The H-resonance lies in the 
neutrino channel for normal mass hierarchy and in the antineutrino channel for the inverted 
hierarchy. The L-resonance lies in the neutrino channel for both the hierarchies [17]. 

The propagation through the resonance L is always adiabatic for the LMA solar param- 
eters. In the case of the H resonance, the "jump" probability depends on the ^13 angle as 
[17]: 



Ph oc exp 



-const sin ^13 



2 / ^m3]_ 



E 



(5) 



Considering Pge = P(fe t'e) and Pee = P(i^e ^e) the survival probabilities, the 
neutrino fluxes arriving at Earth {(t>v) can be written in terms of the fluxes in absence of 
oscillations {(f)°) as: 



= <Pl^Pee + 4>lSl-Pee} 

= rvPee + (t^W^-Pee) 

+ (/.,^ =0°^(1-Pee)+C(l + Pe. 

+ =05 (l-Pee) + 0gJl + ^'e 



(6) 



With ct>i^ = r,, = = <t>lx = 4>i, = 

Using a standard parametrization of the mixing matrix and considering that the neutrino 
bursts do not cross the Earth, the survival probabilities can be expressed in terms of the 
"jump" probability in the H resonance Ph- 



Pee = Ph sin2 cos2 + (1 - Ph) sin2 



Pee = COS 9i2 COS ^13 

for normal hierarchy (Am|2 > 0) and 



(7) 



Pee = sin^ 9i2 COS^ 6*13 

Pee = Ph COS^ 012 COS^ ^13 + (1 - Pff ) sin^ 013 



la ^^o2/i , /I D„^„^„2/) (8) 
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for inverted hierarchy (Amgg < 0). 

The following neutrino oscillations parameters have been used, obtained from other neu- 
trino observations [4, 5]: 

sin2 26'23 = 1. lAmlsl ~ |Am|i| = 3 x IQ-^eF^ 

sin2 6'i2 = 0.3 Am|i = 7 x 10~^eV'^ (9) 

sin^ 013 < 0.02 (at 90% C.L.) Amji > or Am|i < 

According to the adiabaticity conditions on the H resonance and the value of the ^13 
angle, we can consider two limits: the small mixing angle (S) (sin^^ia < 2 x 10~^) and the 
large mixing angle (L) case (sin^ ^^13 > 3 x 10~^). The value of the survival probabilities 
and the neutrino fluxes at Earth in these regions are summarized in Table 2 for both mass 
hierarchies. For intermediate mixing ang le values (2 x lO"'^ < sin^ 6I13 < 3 x IQ-^), the 
"jump" probability Ph depends on ^13 and on the neutrino energy (see e.g. Ref. [11]). 

We can distinguish four extreme cases depending only on the type of mass hierarchy and 
the value of the ^13 mixing angle: 

1. n.h.-L for normal mass hierarchy and large ^13, 

2. n.h.-S for normal hierarchy and small ^13, 

3. i.h.-L for inverted mass hierarchy and large ^13, 

4. i.h.-S for inverted mass hierarchy and small Oi^,. 

For instance, we select the mixing angle sin^ ^13 = 10"^ for the large ^13 case and sin^ ^13 = 
10""^ for the small ^13 case. 

These results have been discussed in Ref. [11]. We recall here the salient feature that 
for normal hierarchy and large mixing angle, the flux on Earth is fully dominated by 

neutrinos in the core of the supernova due to the phenomenon of adiabatic conversion. 
Similarly, for antineutrinos for inverted hierarchy and large mixing angle. 



Limit oscillation 
scenarios 


Survival probabilities 

Pee = P(t'e ^ t'e) 


Neutrino fluxes at Earth 
(l)^=w/o osc, ^cK=with osc 


n.h.-L 

sin^ 6*13 > 3 X 10"^ 


Pee = 0. 
Pee = 0.7 


(t^l = 0.7 (/.g^ + 0.3 ^u.^ 
4>u,+ur+u,+ur = + 0.3 <f>l + 2.7 


n.h.-S or i.h.-S 

sin^ 6'i3 < 2 X 10"^ 


Pee = 0.3 
Pee = 0.7 


= 0.3 + 0.7 r.^ 

= 0.7 + 0.3 

<t}u,+ur+p,+p^ = 0.7 <Pl^ + 0.3 C + 3 


i.h.-L 

sin^ 013 > 3 X 10-^ 


Pee = 0.3 
Pee = 0. 


0., = 0.3 + 0.7 

K+..+l+p^ = 0.7 + C + 2.3 C 



Table 2: Survival probabilities for and i'g and neutrino fluxes at Earth for normal and 
inverted mass hierarchies and large (L) or small (S) mixing angle ^13 in terms of the fluxes 
{(f)°) in absence of oscillations. 
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3 Supernova neutrino signal in LAr 



In a liquid argon TPC supernova neutrinos {E^, < 100 MeV) can be detected through four 
different channels ([11]): 

1. charged-current (CC) interactions on argon 

i^e ^'^Ar -^e- +A' + nN (10) 

i?e^Ar^e+ +A' + nN (11) 

where A' is the remnant nucleus and nN represent emitted nucleons or other debris 
(a's, ...) (if any). At low energy (< 30 MeV) or at low Q^, the dominant reactions are 
Ue '^^Ar —>■ e~ ^^K* and i>e ^^Ar e"*" ^*'C/*. The neutrino energy thresholds of these 
reactions are 1.5 MeV and 7.48 MeV, respectively. The associated photons emitted in 
the de-excitation of the final state nucleus can be used to separate the two reactions. 

2. neutral-current (NC) interactions on argon 

y 40^r +A' + nAr (12) 

where A' is the remnant nucleus and nN represent emitted nucleons or other debris 
(a's, ...) (if any). At low energy (< 30 MeV) or at low Q^, the dominant reaction is 
^^Ar ^^i/' ^^Ar*. The energy threshold of this reaction is 1.46 MeV. The associated 
photons emitted in the de-excitation of the final state nucleus can be used to identify 
the reaction. 

3. elastic scattering (ELAS) on electrons 

(-) - (-) _ 

All these channels have been described in more details in reference [11]. Wc use the 
nuclear cross sections of the CC and NC processes on argon that have been calculated by 
Random Phase Approximation (RPA) for neutrino energies up to 100 MeV including all the 
possible multipoles ([18]). These processes are very important on argon because their high 
cross sections. As already mentioned, it is possible to separate the different channels by 
measuring the associated photons from the K, CI or Ar de-excitation, or by the absence 
of them in the case of elastic scattering. We recall here that the energy dependence of the 
cross-section is oc E"^ for the nuclear processes and oc E for the scattering off electrons. 

Table 3 shows the expected number of neutrino events from a supernova at 10 kpc in a 
100 kton detector. The four detection channels are considered independently and oscillation 
and non-oscillation cases are computed for normal and inverted hierarchies. The supernova 
parameters correspond to scenario I. No threshold on the electron energy has been applied 
since the burst nature of the supernova signal allows to set trigger and detection thresholds 
at the level of the MeV. 
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Scenario I: expected events in 100 kton detector 

= 11 MeV, = 16 MeV, {E,,) = {E^J = 25 MeV 

and luminosity equipartition 



Reaction 


Without 
oscillation 


Oscillation (n.h.) 
Large ^13 Small ^13 


Oscillation (i.h.) 
Large ^13 Small ^13 


ELAS 


1330 


1330 1330 


1330 1330 


l/eCC 
PeCC 


6240 
540 


31320 23820 
1110 1110 


23820 23820 
2420 1110 


NC 


30440 


30440 30440 


30440 30440 


TOTAL 


38550 


64200 56700 


58010 56700 



Table 3: Expected neutrino events in a 100 kton detector from a supernova at a distance of 
10 kpc. 

3.1 Chcirged current events 

In order to understand the dependencies of the expected number of events with the supernova 
parameters and the 613 mixing angle, we perform a scan of these parameters for the UeCC, 
and UgCC processes on argon. While we scan two parameters, the rest are fixed according to 
the values of scenario L 

Figure 1 shows the expected number of UeCC events in a 3 kton detector as a function of 
sin^^i3 and the {E^^J and {E^^) average energies for a normal mass hierarchy. The range of 
the parameters is determined by possible physical regions discussed in Section 2.1, namely 
{Et,^)= (7-18) MeV and {E^J= (15-35) MeV. We can see from these figures that large 
and small mixing angle cases are not symmetric and the results depend on the value of the 
neutrino average energy. fgCC events are very sensitive to the (E^,^) and (E^,^) energies while 
they are of course independent on the (Ep^) value. 

For large ^13 mixing angle, maximal flavor conversion occurs and all v^^r neutrinos convert 
into fe's (see Table 2). Therefore, the number of expected events will only depend on the 
original average energy of Vx neutrinos and not on {Ey^), as can be seen in the figures. For 
small mixing angle the expected flux of u^s at the detector is ~ 30% of original Uf, and 
« 70% of Ux neutrinos produced in the supernova. In this case, the number of fgCC increases 
linearly with {E^,^) and {Ei,^). 

Figure 2 shows the number of Ug^CC events expected in a 3 kton detector as a function of 
sin^ ^13 and the {Ep^ ) and {E^^ ) energies for the case of inverted hierarchy. The H resonance 
lies in the antineutrino channel for the inverted hierarchy and therefore the behavior of the 
i^eCC rates for i.h. is quite similar to the fgCC for n.h. In the first case the variations between 
small and large mixing angle are more pronounced than for fgCC but the event rates statistics 
are smaller. 

The two channels i^gCC and z/gCC are sensitive to the ^13 mixing angle and to the mass 
hierarchy but in different ways. This will help to distinguish the oscillation scenarios. 

3.2 Neutral current events 

Similarly, we perform a scan of the parameters for the NC processes on argon. While we scan 
two parameters, the rest are fixed according to the values of scenario I. 
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Figure 1: Normal mass hierarchy: Variation of the expected number of z^gCC events in a 
3 kton detector as a function of sin^ ^13 and the original supernova neutrino average energies 

{E,:) and {E,^). 




Figure 2: Inverted mass hierarchy: Variation of the expected number of fgCC events in a 
3 kton detector as a function of sin^ ^13 and the original supernova neutrino average energies 
{E,^) and {E,^). 
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Figure 3 shows the variation of the expected events from the NC process with several 
supernova parameters normahzed to a 3 kton detector. 

In the first plot we see the dependence of the NC events with the original average energies 
{Ei,^) and {E^^^) within the ranges defined in Section 2.1, namely {Ei,^)= (7-18) MeV and 
{Ei,^)= (15-35) MeV. Since neutral currents are by definition flavor independent, they are 
not sensitive to the ^13 mixing angle. The number of events increases quickly with average 
energies because of the quadratic dependence of the NC cross section on argon. The variations 
within the range of considered energies of Vx are large while within the range of the (E^,^), 
the variation is much smaller. 

In the second and third plots we scan the average energies as a function of the ratio 
between the electron and non-electron neutrino luminosities [Lg/Lx). As Lg/Lx increases, 
the flux of fe's is increasing with respect to the i^^'s- Accordingly, the number of NC events 
decreases because t'e's are assumed to be less energetic than Vx^s. We can see that the number 

of events increases with {E^^) and {E^^). 

The last plot shows the variations with {Ei,^) and the total binding energy Eb- For low 
(Ei^^) energy, the number of events is almost constant. But as soon as the Vx energy increases, 
the changes with Eb are much bigger, given the quadratic dependence of the cross-section 
on energy. 

All these figures show that the NC events are very sensitive to the different supernova 
parameters independently of the neutrino oscillations. This process hence constitutes an 
excellent probe for the supernova properties independent of the neutrino oscillation physics. 
The combination of the information extracted from other interactions on argon with NC will 
allow to decouple supernova from oscillation physics. 

Elastic events off electrons are also in part sensitive to all neutrino and antineutrino 
flavors. We note however that the rate of NC events is much larger than that of elastic events 
(see Table 3), hence, the statistical power of the nuclear NC events is much larger than that 
of elastic events on electrons. 

3.3 Energy distribution of the charged current and elastic events 

We have presented how the expected rates of the different processes change depending on 
^13, the mass hierarchy and the astrophysical parameters. But also the energy spectra of the 
neutrinos is affected by the value of these parameters. See our reference [11] for examples of 
the variations of the spectra for different values of the ^13 parameter and mass hierarchies. 

Figures 4 and 5 show the variation of the neutrino energy spectra for the elastic scattering, 
fgCC and z/gCC processes as a function of the (E^,^) average energy in the case of large 
mixing angle (sin^ ^13 = 10^^) and normal or inverted mass hierarchies, respectively. The 
total number of events corresponds to the expectations for a 3 kton detector. We have taken 
energy bins of 1 MeV. 

As we consider a large ^13 mixing angle, the fg flux arriving at Earth is ~ 
The energy distribution of the v^CC events depends strongly on the [Ey^) parameter. For the 
inverted mass hierarchy, the main variations as a function of {E^^) are for the f-gCC events. 
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Figure 3: Variation of the expected number of NC events in a 3 kton detector as a function 
of the original supernova U(. and average energies, the ratio between the electron and 
non-electron neutrino luminosities (Le/Lx) and the total binding energy {Eb)- 
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Figure 4: Large mixing angle (sin^ ^13 = 10^) and normal mass hierarchy: Neutrino 
energy spectra for the elastic scattering, UeCC and v'eCC interaction processes on argon as a 
function of the (Ei,^) average energy. The distributions are normahzed to a 3 kton detector. 
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Figure 5: Large mixing angle (sin^ 6*13 = 10^) and inverted mass hierarchy: Neutrino 
energy spectra for the elastic scattering, UeCC and v'eCC interaction processes on argon as a 
function of the (Ei,^) average energy. The distributions are normalized to a 3 kton detector. 
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4 Analysis procedure 



We will consider three different situations. In all cases, we exploit the information from the 
four neutrino detection channels (ELAS, UeCC, i>eCC and NC). When relevant, we compare 
the results with the case of not having access to NC events, which is experimentally the most 
challenging signature. 

1. Study of the different oscillation scenarios without knowledge of the su- 
pernova parameters: We investigate the sensitivity to the ^13 and mass hierarchy 
parameters without any assumption on the supernova parameters. Quantitative results 
on the determination of the sin^ ^13 mixing angle are obtained. This is described in 
Section 5. 

2. Study of the supernova parameters assuming that the oscillation parame- 
ters are known: Assuming that we know the value of the oscillation parameters from 
terrestrial experiments, we study the possible determination of the five supernova pa- 
rameters Eb, {E,yJ, {EpJ, {Ey^) and L^/Lx- Two cases are assumed: (a) the ^13 is 
large (sin^ ^13 = 10~^) and has been measured with a 10% precision from long-baseline 
experiments; (b) the ^13 is small (sin^ ^13 < 10~^) and an upper limit on its value has 
been determined. This is described in Section 6. 

3. Study of the supernova parameters without any knowledge on the neutrino 
oscillation parameters: We investigate the information that can be extracted from 
supernova neutrinos if a supernova explosion occurred nowadays when no information 
about the oscillation parameters ^13 and mass hierarchy (sign[Am32]) is available. This 
is described in Section 7. 

As was shown in the previous section, the simultaneous observation of ELAS, fgCC, v^GC 
and NC events on argon enables us to extract information about the oscillation parameters 
^13 and the type of mass hierarchy (sign[Am|2]) and the five supernova parameters. 

In order to determine all these parameters quantitatively we use a method. This 
analysis procedure includes all the uncertainties and parameter correlations. Wc define a set 
of reference values for the relevant supernova and oscillation parameters that we consider as 
"true". Different values of the ^13 angle are chosen in the analysis for normal and inverted 
mass hierarchies. Two sets of reference values of the astrophysical parameters are studied 
according to scenarios I and II described in Table 1. 

To see how well we can discriminate every "true" scenario, we minimize the in the space 
of the five astrophysical parameters and the two oscillation parameters. The x^ function is 

defined as: 



where A^bms is the number of bins considered in every histogram; Ni is the number of events 
per bin i; a is the statistical error associated to the bin %\ x corresponds to the set of 
parameters to be determined x = {Eb, {Ey^), {EpJ, (E^^), L^/Lx, sin^^is, sign[Ar7T,|2]}, 
which are freely varying and x^ are the "true" values of these parameters. When computing 
the x^ for one parameter, the others are left free. In this way, we can naturally take into 




(14) 
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account the uncertainty on the given parameter introduced by the lack of knowledge on the 
others. This is particularly relevant in the discussion of the decoupling of supernova and 
oscillation parameters. 

The total is computed as the sum of the minimized corresponding to every detection 
channel: 

2 _ 2 I 2 I 2 I 2 /I c-\ 

Xtotal — XeLAS + Xu^CC + XPeCC + XnC U^j 

The fit is performed with the MINUIT [19] package, and is expected to get back the 
same values of the parameters, starting from the reference distributions. At each iteration, a 
different set of parameters is probed, and with the same procedure used to get the reference 
histograms. The condition of the hierarchy between the average energies, confirmed by various 
supernova simulations, is enforced during the minimization: 

(Ke) < (^Pe) < {EuJ (16) 

Otherwise, the ranges of the parameters are left free. 

We compare the energy distributions of the elastic scattering, v^CC and v'gCC neutrino 
events from the "true" scenarios with those obtained with all possible values of the param- 
eters, taking into account the oscillations inside the SN matter. A bin width of 1 MeV is 
considered in the spectra. For the NC channel, the energy distribution cannot be determined 
due to the neutrino presence in the final state. We only consider the number of NC events 
in the fit. 

In order to compute the precision of the determination of the parameters, we consider 
two methods: (1) a onc-dimcnsional "scan" of a given parameter; the other variables are left 
free and minimized at each step; the rcsp. 1, 2, 3 sigmas are given by resp. Xmin + 1' +^ 
+9. (2) a two-dimensional "scan" of a two-parameter plane; the other variables are left free 
and minimized at each point in the plane; the resp. 68%, 90%, 99% C.L. are given by resp. 
xLn + 2-4, +4.6 and +9.2. 

5 Study of the different oscillation scenarios without knowl- 
edge of the supernova parameters 

In this first section we investigate the possible oscillation scenarios without any assumption 

on the supernova properties. 

As was shown in table 2, the value of the survival probabilities Fee and Pee is constant in 
two extreme regions of the ^13 parameter: the large mixing angle case (sin^ ^13 > 3 x 10~^) 
and the small mixing angle case (sin^^ia < 2 X IQ-^). It means that the neutrino fluxes do 
not change with ^13 in these regions. On the other hand, if the angle is in an intermediate 
range 2 x 10~^ < sin^ ^13 < 3 x 10~^, the fluxes depend on the ^13 value and on the neutrino 
energy. Hence, different energy spectra and event rates are expected for every 6*13 value. 

For this reason we first analyze the true small and large mixing angle scenarios and study 
the possibility to distinguish the mass hierarchy and put a bound on the ^13 value. Then, 
we investigate different values of the angle in the intermediate region and we determine the 
statistical accuracy with which it can be measured. 
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5.1 Small and large mixing angle extreme cases 

We consider the four extreme oscillation cases (see Section 2.2) as possible "true" values 
and we use the information of the neutrino detection channels to discriminate between them. 
Because we assume no knowledge on the original supernova neutrino parameters, we let them 
free while performing the minimization. For example, for every value of the ^13 and a 
fixed hierarchy, we vary the supernova parameters in order to obtain the set of values that 
minimize the function for this given value of the mixing angle and mass hierarchy. The 
value result of this fit can then be studied as a function of ^13 and the mass hierarchy. 

A summary of the results on the determination on the ^13 mixing angle for the different 
oscillation cases can be found in Table 4. We include the results for a 3 kton and a 100 kton 
detector and supernova scenarios I and II. 



Determination of ^13 and mass hierarchy (SN parameters free) 


"True" 
osc. 
case 


3 kton LAr 
sin^ ^13 and assumed hierarchy 

SN seen I SN seen II 


100 kton LAr 
sin^ 9i3 and assumed hierarchy 
SN seen I SN seen II 


Events 
type 


n.h.-L 


> 2.1x10^"* n.h. 
excluded i.h. 


> 6.2x10"" n.h. 
excluded i.h. 


> 4.5x10""* n.h. 
excluded i.h. 


> 1.8x10"* n.h. 
excluded i.h. 


all 
all 


> 7.5x10"" n.h. 
excluded i.h. 




> 2.5x10"* n.h. 
excluded i.h. 


> 9.4x10"" n.h. 
excluded i.h. 


no NC 

no NC 


n.h.-S 


< 2.4x10"^ n.h. 

< 3.8x10"^ i.h. 




< 4.3x10"*^ n.h. 

< 5.9x10"^ i.h. 


< 1.2x10"" n.h. 

< 4.3x10"*^ i.h. 


all 
all 


< 3.7x10"'' n.h. 

< 4.1 xlO--"^ i.h. 




< 7.3x10-" n.h. 

< 7.3x10-^' i.h. 


< 3.2x10-" n.h. 

< 1.1x10-'' i.h. 


no NC 
no NC 


i.h.-L 


exchided n.h. 
> 2.3x10-^^ i.h. 




excluded n.h. 
> 5.0x10-'' i.h. 


excluded n.h. 
> 2.3x10-" i.h. 


all 
all 


< 2.8x10"" n.h. 




excluded n.h. 
> 4.9x10-'' i.h. 


< 2.0x10"" n.h. 


no NC 
no NC 


i.h.-S 


< 2.4x10"" n.h. 

< 3.8x10"'^ i.h. 




< 5.9x10"** n.h. 

< 4.0x10"^ i.h. 


< 1.2x10"" n.h. 

< 4.1x10"^ i.h. 


all 

all 


< 3.7x10-'' n.h, 

< 4.()X lO-" i.h. 




< 7.2x10-" n.h. 

< 7.4x10-" i.h. 


< 3.2x10-" n.h 

< l.lxl()-J- i.h. 


no NC 
no XC 



Table 4: Estimated limit on the ^13 mixing angle at la for different "true" oscillation cases 
with a 3 and a 100 kton detectors. We compare the results considering as reference values 
for the supernova neutrino parameters scenarios I (hierarchical) and II (non-hierarchical). 



5.1.1 Supernova scenario I and 3 kton detector 

We first discuss results for a 3 kton detector. Plots in figure 6 show the variation of the 
minimized x^ value with the sin^ ^13 parameter and the mass hierarchy. In the top right 
corner of the figure is indicated in bold the "true" scenario considered (the x° parameters). 
The left part of a given plot corresponds to the results of the fit under the assumption of 
normal hierarchy and the right part for an inverted hierarchy. 

Solid lines are computed considering the four supernova neutrino detection channels. 
Dashed lines do not take into account the NC processes, only ELAS and CC events. The 
one- and two- sigma levels are shown by horizontal lines. The four plots on resp. the top 
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Figure 6: Extreme mixing angle cases and 3 kton detector: value of the fit as 
a function of sin^ ^13 for a 3 kton detector. Curves on the left (right) part of every plot 
correspond to the results of the fit under the assumption of normal (inverted) hierarchy. Four 
limiting oscillation scenarios (n.h.-L, n.h.-S, i.h.-L, i.h.-S) are considered as "true" and the 
reference values of the supernova parameters correspond to scenario I (top) and II (bottom). 
Solid curves are computed using the information of the four supernova neutrino detection 
processes and the dashed line does not take into account NC events. 
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(bottom) are obtained using as true supernova parameters those from resp. scenario I (II). 
The shaded region illustrates the excluded region by the results of reactor experiments [5] . 

We see that normal and inverted hierarchies are indistinguishable for small "true" ^13 
mixing angle. The results of the fit are similar and only an upper bound on sin Oi^ can 
be set. This limit is sin^ ^13 < 2.4 x 10~^ for assumed n.h. and sin^ ^13 < 3.8 x 10~^ for 
assumed i.h. 

If the "true" mixing angle is large, we are able to distinguish among hierarchies. Assuming 
the "true" hierarchy is normal and the value of ^13 is large (n.h.-L), we could put a lower limit 
on the angle (sin^ ^13 > 2.1 x 10~^) for an assumed n.h. The assumed i.h. is excluded. 
The same is possible for the case of inverted hierarchy being the limit sin^ ^13 > 2.3 x 10~^. 
The assumed n.h. is excluded. 

We illustrate the importance of the neutral current events. If we do not include the 
information given by NC events, we loose sensitivity in the large mixing angle cases. The 
limit in the case of n.h.-L can be set only at the level of la and no limit can be determined 
for i.h.-L. Moreover, without NC events, i.h.-L could be misidentified as n.h.-S. 

5.1.2 Supernova scenario II and 100 kton detector 

If we consider the supernova scenario II, we see that the results of the fit obtained with a 
3 kton detector are quite marginal. In this case the average energies of Ux are very close 
to those and i^e- It is impossible to distinguish among hierarchies and no bounds on the 
mixing angle can be set. Only a lower limit is obtained for the case of n.h.-L (sin^ ^13 > 6.2 
x 10-5). 

In order to be sensitive to the oscillation parameters even in the scenario II we consider 
a 100 kton detector (see Figure 7). In this case the sensitivity to the parameters is almost 
the same for scenarios I and II. We can appreciate that the neutral current events are crucial 
even with such a big detector. They are specially important to discriminate the large mixing 
angle regions if the average energies of the neutrinos are close, as we already discussed. 

5.2 Intermediate value of the mixing angle 

If the ^13 mixing angle is in the intermediate range (2 X 10-^ < sin^ 6'i3 < 3 X IQ-^), maximal 
sensitivity to the angle is achieved and measurements of the value are possible in this region, 
since one expects an energy dependent effect in the relevant energy region of the supernova 
neutrinos (sec e.g. Ref. [11]). 

If we consider a mixing angle sin^ ^13 = 10~^, the fit gives the rcsTilts shown in figure 
8. True normal (n.h.-i) and inverted (i.h.-i) hierarchies for true sin^ 613 = 10"^ and the 
two supernova scenarios are studied. Solid lines correspond to the contribution of the four 
neutrino detection channels while dashed lines only consider elastic and CC events. The one- 

and two- sigma levels are shown by horizontal lines. 

A determination of the angle is possible at la level for scenario I with a 3 kton detector 
considering all neutrino processes. We obtain the following la intervals and mass hierarchies 
for a 3 kton detector: 

5.8 X 10"^ < sin^ 6*13 < 1.9 x lO""^ {assumed n.h.) (17) 
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Figure 7: Extreme mixing angle cases and 100 kton detector: value of the fit as a 

function of sin^ ^13 for a 100 kton detector. Curves on the left (right) part of every plot cor- 
respond to the results of the fit under the assumption of a normal (inverted) hierarchy. Four 
limiting oscillation scenarios (n.h.-L, n.h.-S, i.h.-L, i.h.-S) are considered as "true" and the 
reference values of the supernova parameters correspond to scenario I (top) and II (bottom). 
Solid curves are computed using the information of the four supernova neutrino detection 
processes and the dashed line does not take into account NC events. 
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Figure 8: Intermediate mixing angle region (sin^ ^13 = 10~^): value of the fit as 
a function of sin^ ^13 for a 3 kton detector. Curves on the left (right) part of every plot 
correspond to results of the fit for assumed normal (inverted) hierarchy. Two true oscillation 
scenarios are considered (n.h. and i.h.) and the true values of the supernova parameters 
correspond to scenario I (left) and II (right). Solid curves are computed using the information 
of the four supernova neutrino detection processes and the dashed line does not take into 
account NC events. 
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3.3 X 10~^ < sin^ < 1.8 x 10~^ {assumed i.h.) (18) 

The wrong hierarchies are excluded in both cases. 

However, the statistics achieved with such a mass are not enough to put any constraint 
on the angle if the supernova parameters correspond to scenario II. Figure 9 compares 
the determination of the sin^ ^13 parameter using a 3 and a 100 kton detector. The of the 
fit is plotted as a function of the sin^ ^13 value. 
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Figure 9: Determination of the sin^ ^13 parameter in the intermediate region. We have 
considered as "true" values sin^ 61^ = 10~^, normal (left) and inverted (right) hierarchies and 
the supernova parameters of scenario I. The solid curve corresponds to the results of the 
fit obtained with a 100 kton detector and dashed line shows the contribution of a 3 kton. 
The one- and two- sigma levels are shown by horizontal lines. 

If we consider a 100 kton detector, much more statistically accurate measurements can 
be performed at la and 2a levels: 



9.1(8.1) X 10-5 < sin2 6'i3 < 1.1(1.3) x 10" 



8.9(7.4) X 10"^ < sin^0i3 < 1.1(1.3) x 10 
The wrong hierarchies are excluded in both cases. 



-4 



la{2a) {assumed n.h.) (19) 
la{2a) {assumed i.h.) (20) 



The precision on the measurement of the mixing angle will depend on the exact value 
of 9is. Considering different reference values in the intermediate region, Figure 10 shows 
the superposition of different determinations of the sin^ 9is parameter, for the cases of true 
normal and inverted hierarchies normalized to a 3 kton detector. The angle can be measured 

at la in both cases. 

The estimated limits and ranges of sin^ ^13 are summarized in table 5 for true normal 
and inverted hierarchies. If the sin^ ^13 angle is between (4x10"^ - 1x10"^) (n.h.) or 
(7xl0~5 - 3xl0~^) (i.h.), it could be constrained in a certain range. Otherwise, upper or 
lower limits on its value can be set. 
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Figure 10: Determination of the sin^ ^13 parameter in the intermediate region with a 3 kton 
detector. We have considered different "true" values for sin^ ^13 and normal and inverted 
hierarchies. The reference supernova parameters correspond to scenario I. The one- and two- 
sigma levels are shown by horizontal lines. 



3 kton detector 


"True" value 

sin^ 6I13 


Limit or range of sin'^ ^13 
at Icr level 


True n.h. 


n.h. assumed 


1. X 10-^ 


sin^6li3 < 3.4 x IQ-^ 


2. X 10^5 


sin^6'i3 < 4.6 x IQ-^ 


4. X 10-^ 


1.6 X 10-5 < sin^6li3 < 7.2 x 10-5 


6. X 10^5 


2.8 X 10-5 < sin^6li3 < 9.4 x 10-5 


8. X 10-^ 


4.2 X 10-5 < sin''6li3 < 1.3 x 10-* 


1. X 10-'' 


5.8 X 10-5 < sin^6li3 < 1.9 x 10-* 


2. X 10--^ 


sin^6'i3 > 1.1 X IQ-'' 


3. X 10-'' 


sin^6ii3 > 1.5 X 10-4 


True i.h. 


i.h. assumed 


2.5 X 10-*^ 


sin^6ii3 < 4.0 X 10-5 


2.5 X 10-5 


sin^6ii3 < 6.8 X 10-5 


5.0 X 10-5 


sin^ 9i3 < 9.9 X 10-5 


7.7 X 10-5 


1.3 X 10-5 < sin^6li3 < 1.4 x 10-* 


1.0 X lO-'' 


3.3 X 10-5 < sin'^6li3 < 1.8 x 10-* 


1.3 X lO-'' 


5.4 X 10-5 < sin^ 613 < 2.2 x 10-* 


1.7 X 10-'^ 


8.5 X 10-5 < sin^6li3 < 3.4 x 10-* 


2.0 X 10--* 


1.1 X 10-* < sin^^^is < 5.3 X 10-* 


3.3 X 10--* 


sin^6li3 > 1.7 X 10-* 



Table 5: Estimated limit on the ^13 mixing angle at la level for different "true" values of the 
angle and mass hierarchies in the intermediate range. 
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6 Study of the supernova parameters assuming that the os- 
cillation parameters are known 



In this section we study the information about the astrophysical parameters that can be 
extracted from the detection of supernova neutrinos assuming that the oscillation parameters 
sin^ 013 and the mass hierarchy have been determined from terrestrial experiments, essentially 
from long-baseline experiments. 

Considering the range of sensitivity of ^13 of future accelerator experiments (see e.g. [20]), 
we can study two cases: 

1. If ^13 angle is measured by long-baseline experiments: If the true value of the 
013 angle is large (sin^ ^13 > 10^^) then future superbeams should be able to measure 
it (see e.g. [20]). We assume that the angle will be determined with a precision of 10%. 

Therefore, we study how well we can constrain the supernova parameters under the 
external constraint sin^ ^13 = 10~^ ± 10~^, considering that the mass hierarchy is also 
known. 

2. If an upper bound on Bis angle is set by long-baseline experiments: If the 

angle is very small, then future long-baseline experiments might not have the possibility 
to measure this angle and will only place an upper limit on the value. As a second 
example, we consider that the limit sin^ ^13 < 10~^ has been set and we investigate the 
consequences of this limit into the determination of the supernova parameters. 

6.1 If ^13 angle is measured by long-baseline experiments 

We assume that the value of the true 613 mixing angle has been determined by long-baseline 
experiments to be sin^ ^13 = 10^'^ it 10^^ and we study the cases of true normal and inverted 
hierarchies. We will compare the precisions in the determinations of the supernova parameters 
with a 3 and 100 kton detectors. In the fits, we let free the supernova parameters {Eb, {E^J, 
(Ep^), {EyJ), Le/Lx) and we perform a minimization in order to obtain the allowed regions 
in different ID parameter regions and 2D parameter planes. The expected accuracies at 90% 
C.L. to be achieved with a 3 and 100 kton detectors are summarized in Table 6. 

Figure 11 shows on the top the value of the fit as a function of the supernova param- 
eters for a 3 kton detector for true normal hierarchy. Figures 12 correspond to the inverted 
hierarchy case. The bottom plots arc the 68%, 90% and 99% C.L. allowed regions for the 
astrophysical parameters with a 100 kton detector. The crosses on the plots indicate the 
values of the parameters taken as reference for the fit, which correspond to scenario I. 

The plots Le/Lx vs (E^) show the luminosity of every flavor as a function of its average 
energy. We see that with the statistics given by the 3 kton detector is not possible put any 
constraint on the {E^^) and (-Ep^) energies or on Le/Lx- However, the (E^^) energy can be 
determined very precisely. The accuracy that we can achieve at 90% C.L. is /S.{Ei,^) / [Ei,^)^ 
3%. 

The plot Eb vs L^/Lx shows the strong correlation between these two variables. Using 
a 100 kton detector, all the variables can be determined in a certain range. The best results 
are obtained for [Ei,^ ) because we are considering the large ^13 mixing angle. In this case the 
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With constraint sin^ 


9^■^ = 10~^ 


lb 10 ^ from terrestrial experiment 




ueuecror mass 


True hierarchy 


SIM cjppn 






(E,J 


AEg 

Eb 




3 kton 


n.h. 


I 






~ 3% 


~ 23% 




3 kton 


i.h. 


I 


~ 46% 




~ 7% 


~ 17% 


~ 60% 


3 kton 


n.h. 


II 




~ 26% 


~ 4% 


~ 28% 




3 kton 


i.h. 


II 




~ 36% 


19% 


~ 20% 




100 kton 


n.h. 


I 


~ 14% 


- 4% 


< 1% 


~ 2% 


~ 11% 


100 kton 


i.h. 


I 


~ 5% 


- 9% 


~ 1% 


~ 2% 


~ 9% 


100 kton 


n.h. 


II 


~ 9% 


- 3% 


< 1% 


~ 4% 


~ 12% 


100 kton 


i.h. 


II 


~ 6% 


~ 6% 


~ 1% 


~ 2% 


~ 32% 



Table 6: Expected accuracies at 90% C.L. in the determination of the supernova parameters 
using the neutrinos measured with a 3 kton and a 100 kton detector. We have assumed that 
the mass hierarchy is known and the ^13 mixing angle has been measured by long-baseline 
experiments, being equal to sin^ ^13 = 10~^ it 10~^. Supernova scenarios I and II are tested. 



I'e's have oscillated into i/^/s and i^t-'s (so J^^'s) and vice versa. Therefore, the most important 
reaction on argon (z/gCC) will be mainly sensitive to the original {E^^). 

In the inverted hierarchy fraction of the flux arrives at the detector hence in 

this case we are sensitive to the {E^J energy and a bit less to (E^^), but both variables can 
be estimated with accuracies of ~ 46% and ~ 7%, respectively, for a 3 kton detector. 

If we consider true values those given by scenario II (see table 1), the corresponding results 
for the fit are presented in figures 13 (n.h.) and 14 (i.h.). For the normal hierarchy case we 
sec that results are quite similar to those obtained for scenario I. Only the (E,^^) parameter 
can be measured with a 3 kton detector. However, the inverted hierarchy case shows a worse 
situation. The Ue and average energies are closer for scenario II and it makes more difficult 
the determination of a region in the plane {{Ei,^), [En^)). 

For a 100 kton detector, scenarios I and II give quite similar results except for the 
L^/Lx parameter in the case of i.h., which is worse determined for scenario II (~ 32%) 
than for scenario I (~ 9%) (see table 6). 

6.2 If an upper bound on 6'i3 angle is set by long-baseline experiments 

Another possibility is that future long-baseline neutrino experiments will not be sensitive 
enough to measure sin^ ^13 and will place an upper limit on its value. Assuming the limit 
sin^ ^13 < 10~^, we study the possible determination of the supernova parameters. 

In this case the large mixing angle region for both hierarchies is excluded and we only 
have the possibility that the ^13 angle be in the intermediate (n.h.-i or i.h.-i) or in the small 
(n.h.-S or i.h.-S) mixing angle regions. 

Putting the condition of sin^ Oi^, < 10~^ and fixing the mass hierarchy, we perform a 
minimization letting at each step all variables free, included the ^13 angle. The upper bound 
on the 613 angle is inserted as an appropriate term in the function. We study two cases 
for the "true" data: normal and inverted hierarchy. 

Table 7 shows the expected accuracies at 90% C.L. in the determination of the supernova 
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Figure 11: (Top) value of the fit as a function of the supernova parameters for a 3 kton 
detector, assuming that the ^13 mixing angle has been measured with a precision of 10% 
(sin^ 6*13 = 10-3 ± 10-"^) and the mass hierarchy is normal {Am^^ > 0). (Bottom) 68%, 
90% and 99% C.L. allowed regions for the supernova parameters with a 100 kton detector. 
Crosses indicate the value of the parameters for the best fits. 



24 




100 kton LAr 









— 68% C.L 








--- 90% C.L 








99% C.L 












_ + 






[ 










i 16 - 

f " - 
12 -■; 

10 



< E- > (MeV) 




25 30 35 

< E^^ > (MeV) 





— 68% C.L 




--- 90% C.L 




99% C.L 


L, , , , 1 , , , 


!, , , , 1 , , , 



< E„ > (MeV) 



Figure 12: (Top) value of the fit as a function of the supernova parameters for a 3 kton 
detector, assuming that the ^13 mixing angle has been measured with a precision of 10% 
(sin^ 6I13 = 10"^ ± 10-"^) and the mass hierarchy is inverted (Am^^ < 0). (Bottom) 68%, 
90% and 99% C.L. allowed regions for the supernova parameters with a 100 kton detector. 
Crosses indicate the value of the parameters for the best fits. 
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sin2 6'i3 = 10-3 ± 10-4 and n.h. (seen II) 
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Figure 13: (Top) value of the fit as a function of the supernova parameters for a 3 kton 
detector, assuming that the ^13 mixing angle has been measured with a precision of 10% 
(sin^ 013 = 10-3 ± 10-"^) and the mass hierarchy is normal {Arri^i > 0). The reference values 
taken for the neutrino average energies are the ones of scenario II (see table 1). (Bottom) 68%, 
90% and 99% C.L. allowed regions for the supernova parameters with a 100 kton detector. 
Crosses indicate the value of the parameters for the best fits. 
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sin2 6'i3 = 10-3 ± 10-^ and i.h. (seen II) 
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Figure 14: (Top) value of the fit as a function of the supernova parameters for a 3 kton 
detector, assuming that the ^13 mixing angle has been measured with a precision of 10% 
(sin^ (9i3 = 10-3 ± 10"^) and the mass hierarchy is inverted (Aml^ < 0). The reference values 
taken for the neutrino average energies are the ones of scenario II (see table 1). (Bottom) 68%, 
90% and 99% C.L. allowed regions for the supernova parameters with a 100 kton detector. 
Crosses indicate the value of the parameters for the best fits. 
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parameters for 3 and 100 kton detectors and supernova scenarios I and II. 



With constraint sin'^ ^13 < 10 from terrestrial experiment 


Detector mass 


True hierarchy 


SN seen. 






(E.J 


Eb 


MLe/LJ 
{Le/LJ 


3 kton 


n.h. 


I 


~ 30% 


- 30% 


~ 7% 


~ 26% 




3 kton 


i.h. 


I 


~ 52% 


~ 30% 


- 8% 


~ 17% 




3 kton 


n.h. 


II 




- 25% 


~ 16% 


~ 27% 




3 kton 


i.h. 


II 




- 25% 


~ 23% 


~ 17% 




100 kton 


n.h. 


I 


~ 6% 


~ 4% 


< 1% 


~ 1% 


~ 11% 


100 kton 


i.h. 


I 


~ 5% 


~ 4% 


< 1% 


~ 2% 


~ 9% 


100 kton 


n.h. 


II 


- 9% 


~ 5% 


< 1% 


~ 2% 


~ 35% 


100 kton 


i.h. 


II 


~ 8% 


~ 4% 


~ 1% 


~ 2% 


~ 37% 



Table 7: Expected accuracies at 90% C.L. in the determination of the supernova parameters 
using the neutrinos measured with a 3 kton and a 100 kton detector. We have assumed that 
the mass hierarchy is known and an upper limit on the mixing ang le of sin2 6*13 < 10"^ has 
been set by long-baseline neutrino experiments. Supernova scenarios I and II are tested. 



Figures 15 and 16 show on the top the value of the fit as a function of the supernova 
parameters for a 3 kton detector considering that an upper limit on the value of ^13 has been 
set (sin^ ^13 < 10"*^) and for normal or inverted mass hierarchy, respectively. The bottom 
plots are the allowed regions at 68%, 90% and 99% C.L. on different supernova parameters 
planes for a 100 kton detector. A value of sin^ ^13 = 10~^ was taken as reference for the fit, 
i.e. in the small mixing angle region. However, the fitted angle can be in both intermediate 
and small mixing angle regions. We have checked that if the "true"' value of the angle is 
sin^ 013 = 10~^, the results of the fit are only slightly modified. 

In this case the poor determination of the parameters obtained with 3 ktons indicates 
that statistics are not enough for providing good results. However, a 100 kton detector will 
solve these problems. 

Normal and inverted hierarchy cases give different results because they depend on the 
value of the angle. For sin^ ^13 < 10~® it is impossible to distinguish between both hierarchies. 
But if the angle is the intermediate region (2 X 10-6 < sin^^ia < 3 x 10"^) there IS a 
dependence with the neutrino energy and the results depend on the particular value of ^13. 
This is specially important in the n.h. case, where the survival probability changes with 
energy in the neutrino channel, which is the most sensitive for LAr detectors. 

With a 100 kton detector, the expected accuracies are similar for normal and inverted 
hierarchies, as shown in table 7. Thanks to the high statistics available with a very massive 
detector, we can fit very precisely the data and the results are essentially equal for both mass 
hierarchies, as expected in the small ^13 region. 

Considering the parameters of scenario II as reference values, the corresponding 2D al- 
lowed regions are presented in figures 17 and 18, for both mass hierarchies. In this case the 
determination of {Ey^) and {Ey^) is more complicated due to the non-hierarchical scenario. 
Likewise, L^lLx presents a big uncertainty even with a 100 kton detector. 



28 



sin2 6'i3 < 10-^ and n.h. 
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Figure 15: (Top) value of the fit as a function of the supernova parameters for a 3 kton 
detector, assuming that an upper hmit on the value of the ^13 mixing angle has been set 
(sin^ 9i3 < 10"^) and the mass hierarchy is normal {Am^i > 0). (Bottom) 68%, 90% and 
99% C.L. allowed regions for the supernova parameters with a 100 kton detector. Crosses 
indicate the value of the parameters for the best fits. 
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sin^6>i3 < 10"^ and i.h. 
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Figure 16: (Top) value of the fit as a function of the supernova parameters for a 3 kton 
detector, assuming that an upper hmit on the value of the ^13 mixing angle has been set 
(sin^ 6113 < 10"^) and the mass hierarchy is inverted {l^m^i < 0). (Bottom) 68%, 90% and 
99% C.L. allowed regions for the supernova parameters with a 100 kton detector. Crosses 
indicate the value of the parameters for the best fits. 
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sin^6'i3 < 10 ^ and n.h. (seen II) 
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Figure 17: (Top) value of the fit as a function of the supernova parameters for a 3 kton 
detector, assuming that an upper hmit on the value of the ^13 mixing angle has been set 
(sin-^ 013 < 10"^) and the mass hierarchy is normal {Am^^ > 0). The reference values taken 
for the neutrino average energies are the ones of scenario II (see table 1). (Bottom) 68%, 
90% and 99% C.L. allowed regions for the supernova parameters with a 100 kton detector. 
Crosses indicate the value of the parameters for the best fits. 
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sin2 6'i3 < 10-^ and i.h. (seen II) 
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Figure 18: (Top) value of the fit as a function of the supernova parameters for a 3 kton 
detector, assuming that an upper hmit on the value of the ^13 mixing angle has been set 
(sin^ 9i3 < 10"'') and the mass hierarchy is inverted (Am^^ < 0). The reference values taken 
for the neutrino average energies are the ones of scenario II (see table 1). (Bottom) 68%, 
90% and 99% C.L. allowed regions for the supernova parameters with a 100 kton detector. 
Crosses indicate the value of the parameters for the best fits. 
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7 Study of the supernova parameters without any knowledge 
on the neutrino oscillation parameters 



We consider four cases as "true" values of the oscillation parameters: n.h.-L, i.h.-L, n.h.-S, 
i.h.-S. Without assuming any value of the mixing angle and mass hierarchy, we perform a ^ 
minimization letting all parameters {Ebi {Ei/^)-, {Ep^), {E^^), L^/Lx, sin^^is, sign[Am|2]} 
freely vary and we study how well we can determine them in the case of a 100 kton detector. 
For the 3 kton detector, the results can be rescaled accordingly. 

Figures 19 and 20 show the allowed regions at 68%, 90% and 99% C.L. of the oscillation 
and supernova parameters. The reference values of the parameters correspond to scenario I. 

The results show that with large statistics it is possible to decouple the supernova and 
oscillations parameters and determine the parameters with high precision, as presented in 
Table 8. The expected accuracies at 90% C.L. are listed for every astrophysical parameter, 
supernova scenario and "true" oscillation case considered. The total binding energy {Eb) 
has an error 1-4%, the average energy of electron neutrinos {Eu^) at the core is determined 
with an error 5-21%, the average energy of electron antineutrinos (Ep^) at the core, with 
an error 4-9%, the average energy of other (anti)neutrinos {E^J at the core is determined 
with an error ~ 1%, and the relative luminosities of the electron and non-electron flavor 
neutrinos L^/Lx, with an error 9-37%. The errors with a 3 kton detector are roughly 6 times 
larger than those with a 100 kton detector, hence, many of the parameters would be poorly 
determined in this case. 



100 kton detector. All parameters free. 


"True" osc. case 


SN seen. 




A (AW) 




AEb 
Eb 


A(Le/L,) 
(Le/L:.) 


n.h.-L 


I 


~ 17% 


~ 4% 


< 1% 


~ 2% 


~ 11% 


i.h.-L 


I 


~ 5% 


~ 9% 


< 1% 


~ 2% 


~ 9% 


n.h.-S 


I 


~ 6% 


~ 4% 


< 1% 


~ 1% 


~ 11% 


i.h.-S 


I 


~ 5% 


~ 4% 


< 1% 


~ 2% 


^ 9% 


n.h.-L 


II 


~ 21% 


~ 3% 


< 1% 


~ 4% 


~ 14% 


i.h.-L 


II 


~ 6% 


~ 9% 


~ 1% 


~ 2% 


~ 34% 


n.h.-S 


II 


~ 11% 


~ 5% 


< 1% 


~ 2% 


~ 35% 


i.h.-S 


II 


- 8% 


~ 4% 


~ 1% 


~ 2% 


~ 37% 



Table 8: Expected accuracies at 90% C.L. in the determination of the supernova parameters 

using the neutrinos measured with a 100 kton detector. No conditions on the ^13 angle and 
mass hierarchy have been considered. Four oscillation scenarios have been studied as possible 
"true" cases. Supernova scenarios I and II are tested. 

Figure 21 shows the value from the minimization as a function of the supernova 
parameters for a 100 kton detector, considering scenario I (left) and II (right) as "true" 
supernova parameters. We compare the results obtained fixing the mixing angle to sin^ ^13 
= 10~^ ± 10~^ and normal mass hierarchy with the case of leaving the mixing angle free. 
We have chosen the case where the effect of the neutrino mixing is largest, namely when the 
"true" oscillation parameters correspond to normal mass hierarchy and large mixing angle. 
For other oscillation parameters, in particular for small mixing, the effect is invisible. 

We see that the expected accuracies are similar to the ones obtained constraining the 
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n.h.-L and n.h.-S (all parameters free) 
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Figure 19: Correlation between oscillation and supernova parameters: 68%, 90% and 99% 
C.L. allowed regions for the supernova parameters without any assumption on the oscillation 
parameters for a 100 kton detector. The reference values correspond to the normal mass 
hierarchy, large (left) and small (right) ^13 mixing angle cases. Crosses indicate the value of 
the parameters for the best fits. 
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i.h.-L and i.h.-S (all parameters free) 
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Figure 20: Correlation between oscillation and supernova parameters: 68%, 90% and 99% 
C.L. allowed regions for the supernova parameters without any assumption on the oscillation 
parameters for a 100 kton detector. The reference values correspond to the inverted mass 
hierarchy, large (left) and small (right) ^13 mixing angle cases. Crosses indicate the value of 
the parameters for the best fits. 
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mixing angle and the mass hierarchy, except for the {E^^) parameter and "true" oscillation 
case n.h.-L. In this scenario, the knowledge of the ^13 value improves the determination of 
the {Ey^) energy from 17% to 14% for scenario I and from 21% to 9% for scenario II. 

Hence, the observation of different channels which have different sensitivities to the neu- 
trino flavors allows to decouple the supernova and neutrino oscillation physics, effectively 
allowing these two sectors to be studied independently. With the statistics provided by a 
100 kton detector, a single supernova explosion would allow to determine the parameters of 
the supernova cooling phase quite precisely. 
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Figure 21: Determination of the supernova parameters: value of the fit as a function 
of the supernova parameters for a 100 kton detector. We compare the results obtained fixing 
the angle to sin^ ^13 = 10^'^ ± 10~^ and normal mass hierarchy with the case of leaving all 
parameters free. The "true" supernova parameters correspond to scenario I (left) and to 
scenario II (right). The "true" oscillation parameters are normal mass hierarchy and large 
mixing angle. 



8 Conclusions 

In this paper we have investigated the capabilities of liquid argon TPC detectors to study 
and decouple supernova and neutrino oscillation physics in the event of the detection of a 
single core collapse supernova at an assumed distance of 10 kpc. 

The neutrino physics was determined by the oscillation angle ^13 and the mass hierarchy. 
The supernova properties were summarized in five astrophysical parameters: the average 
energies of the neutrinos emitted from the supernova {{E^^), (Ep^), (E^^)), the total binding 
energy {Eb) and the relative luminosities of the electron and non-electron flavor neutrinos 
(Le/Lx). We considered two specific scenarios (I & II) in order to understand the effects of 
a hierarchical versus non-hierarchical distribution of energies of supernova neutrinos, since 
this issue is still debated. For the supernova scenario II, the differences between the average 
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energies of electron and non-electron neutrinos is of the order of 10%. This is an important 
issue since the interplay between hierarchy between neutrino flavors and neutrino mixing is 
expected to be a striking feature of matter enhanced oscillations in the supernova. In the case 
of degenerate neutrinos, the effect of oscillation is more difficult to observe. To contemplate 
two scenarios allowed us to address this in a quantitative way. 

We have considered the four neutrino detection channels available in argon: elastic scat- 
tering on atomic electrons, charged-current interactions on argon (independently sensitive to 
t'e and Ve neutrinos) and neutral-current interactions on argon. 

In the first part of the work, we have studied the sensitivity to the ^13 and mass hierarchy 
parameters under the assumption of the supernova scenario I. For the true large mixing angle 
(> 3 X 10"^), a mass hierarchy identification is possible and lower limits on the 612, value 
are set. We have shown that the detection of the NC process is important in this context. 
For small mixing angle (< 2 x 10~^), it is not possible to distinguish the mass hierarchy but 
upper limits on can be obtained. If the mixing angle is in an intermediate range (2 x 
10~^ < siv? 9iz < 3 X 10~^), measurements of its value are possible. 

For the supernova scenario II, we found with the statistics provided by a 3 kton detector 
the possibility to determine the ^13 mixing angle is limited. More statistics, such as the 
one obtained with a 100 kton detector, will be mandatory to determine quantitatively the 
oscillation parameter. 

In the second part of our study, we have studied the ability to determine the supernova 
parameters considering two different cases for the ^13 angle: (a) the angle has been determined 
by long-baseline experiments with !=a 10% precision, (b) an upper limit on the mixing angle has 
been set by long-baseline experiments. For a 3 (100) kton detector, the supernova parameters 
are measured with such accuracies: the total binding energy {Eb) has an error 20-30% (2- 
4%), the average energy of electron neutrinos {Ey^) at the core is determined with an error 
30-120% (5-14%), the average energy of electron antineutrinos {Ep^) at the core, with an error 
30—80% (3—9%), the average energy of other (anti)neutrinos {Ei,^) at the core is determined 
with an error 3-23% (< 1%), and the relative luminosities of the electron and non-electron 
fiavor neutrinos Lg/L^, with an error > 100% (10-40%). 

In the last part of our study, we have considered the case where both supernova and 
oscillation parameters are free in the minimization. The precision with which the supernova 
parameters can be determined without knowledge on the neutrino mixing parameter is essen- 
tially the same as when information from terrestrial experiments would be available, except 
in the supernova scenario II and in large mixing angle case where the determination of the 
{En^) energy profits largely from the terrestrial knowledge on ^13. 

In conclusion, the possibility provided by liquid Argon TPCs to observe elastic, charged 
and neutral current channels with different sensitivities to the neutrino flavors allows to 
disentangle the supernova and neutrino oscillation physics, allowing these two sectors to be 
effectively decoupled and studied independently. With the statistics provided by a 100 kton 
detector, a single supernova explosion would allow to determine the supernova cooling phase 
quite precisely. This is true even in the case that the energies of the neutrinos of different 
flavors are almost degenerate, as some recent supernova simulations indicate. 
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